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1.1 Introduction 

Why was there a Big Bang? Why is the universe not featureless and bar- 
ren? Why are there fluctuations in the cosmic microwave background? In 
stark contrast to the convincingly answered "what" questions of cosmol- 
ogy (e.g. What is the age of the universe?, What is the geometry of the 
universe?), these "why" questions may instead evoke a sense of disillusion- 
ment. Is it possible that cosmology's "triumphs" — its answers to the "what" 
questions — are frustratingly inadequate, or worse, incomplete? 

However, what if the "why" questions provide tantalizing hints of the 
ultimate origins of the universe? Then instead of crisis, we encounter an 
amazing opportunity — one that might provide answers to the most enig- 
matic question of all: How did the universe begin? 

Inflation [1] is a daring paradigm with the promise to solve many of these 
mysteries. It has entered its third decade of successfully confronting observa- 
tional evidence and emerged as cosmology's theoretical touchstone. Despite 
its many successes, inflation remains unproven. While skeptics must resort 
to increasingly finely tuned attacks [21 [3], inflation's proponents can only 
cite circumstantial evidence in its favor [Ij. However, a conclusive detection 
of a primordial gravitational wave background (GWB) from inflation would 
be "the smoking gun" [5]. No other known cosmological mechanism mimics 
the GWB's imprint on the cosmic microwave background (CMB). 

New technological innovations poise cosmology at the threshold of an 
exhilarating era — one in which future CMB data will winnow down the 
seemingly boundless "zoo" of cosmological models and test the hypothesis 
that an inflationary expansion of the universe took place in its first moments. 

Inflation's unique imprint on CMB polarization has generated consider- 
able attention from US science policy advisors O [71 El E], who have all 
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enthusiastically recommended measuring CMB polarization. The reason for 
this excitement is clear: inflation explains a host of critical cosmological ob- 
servations, and CMB polarization is the most promising, and perhaps only, 
way to glimpse the GWB. 

This chapter describes how the Cosmic Gravitational Wave Background 
induces a specific type of CMB polarization and describes the first experi- 
ment dedicated to testing this most-promising signature of inflation. This 
experiment, the Background Imaging of Cosmic Extragalactic Polarization 
(BICEP) project, has recently embarked on its third observing season. We 
show preliminary data from the BICEP 's first season obtained with a novel 
polarization modulation mechanism called the "Faraday Rotation Modula- 
tor" . Our discussion ends with a description of exciting new technology with 
the potential to probe infiation down to the ultimate cosmological limit. 

1.2 The inflationary universe 

The CMB has historically been the tool to appraise inflationary cosmology. 
This is not surprising since the CMB is the earliest electromagnetic "snap- 
shot" of the universe, a mere 380,000 years after the Big Bang. As such 
it probes the universe in a particularly pristine state — before gravitational 
and electromagnetic processing. Because gravity is the weakest of the four 
fundamental-forces, gravitational radiation (i.e. the GWB) probes much far- 
ther back: to ~ 10"'^^ seconds after the Big Bang {W^tpi in Planck units). 
The GWB encodes the cosmological conditions prevailing at 10^^ GeV en- 
ergy scales. In contrast, the CMB encodes the physical conditions of the uni- 
verse when radiation decoupled from matter at energy scales corresponding 
to 0.3 eV at t ~ lO^^tpi. As experimentalists, we can exploit the primacy of 
the CMB by using the CMB's surface of last scattering as a "film" to "ex- 
pose" the GWB — primordial reverberations in spacetime itself. Doing so 
will provide a "baby picture" of the infant universe; an "ultrasonic" image 
of the embryonic universe! 

1.2.1 Quantum fluctuations in the inflationary universe 

Infiation posits the existence of a new scalar field (the inflaton) and speci- 
fies an action-potential leading to equations of motion. Quantizing the in- 
fiaton field causes the production of perturbations (zero point fiuctuations) 
|10j . While the infiaton's particle counterpart is unknown, its dynamics as 
a quantum field have dramatic observational ramifications JT]. All viable 
cosmological theories predict a spectrum of scalar (or energy density) per- 
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turbations that can then be tested against CMB temperature anisotropy 
measurements. Inflation predicts the spectrum of scalar perturbations, and 
additionally predicts tensor perturbations (i.e. the GWB). Inflation's unique 
prediction is the GWB, which is parameterized by the tensor-to-scalar ra- 
tio, r. An unambiguous detection of r will reveal both the epoch of inflation 
and its energy scale [5]. If, as theorists have speculated [12], inflation is re- 
lated to Grand Unifled Theories (GUT), then a detection of the GWB also 
will probe physics at energy scales one trillion times higher than particle 
accelerators such as the Large Hadron Collider pj|. 

Both energy density fluctuations (scalar perturbations) and gravitational 
radiation (tensor perturbations; the GWB) produce CMB polarization. The 
two types of perturbations are related in all inflation models, since both are 
generated by quantum fluctuations of the same scalar held, the inflaton [13]. 
The relationship between CMB polarization produced by scalars and ten- 
sors will provide a powerful consistency check on inflation when the GWB 
is detected. Similar relations, using recent detections of the scalar pertur- 
bation spectrum's departure from "scale invariance" (primarily using 
CMB temperature anisotropy), have led to claims of "detection" of infla- 
tion, at least in the popular press [16]. Furthermore, NASA's Wilkinson 
Microwave Anisotropy Probe (WMAP) showed an anti-correlation between 
temperature and polarization at large angular scales, providing additional, 
albeit circumstantial, evidence in favor of inflation |17[ 118]. 

The ultimate test of inflation requires a measurement of the tensor power 
spectrum itself, not only the predicted temperature-polarization correlation 
or the properties of the scalar power spectrum. Given a very modest set of 
external, non-inflation-specific parameters, including a simple cosmological 
chronology (specifying that inflation was followed by radiation domination, 
subsequently followed by matter domination), inflationary models can pre- 
cisely predict the spatial correlations imprinted on the polarization of the 
CMB by the GWB, making it truly "the smoking gun." 

1.2.2 The gravitational wave background: shaking up the CMB 

Scalar metric-perturbations have no handedness, and are therefore said to 
be "parity invariant." While the GWB produces both temperature and po- 
larization perturbations, the temperature perturbations are primarily as- 
sociated with the change in potential energy induced by the gravitational 
waves, whereas the tensor-induced polarization perturbations are associated 
with spacetime stress and strain. The parity-violating polarization signature 
exists only if cosmological gravitational waves exist, as flrst demonstrated 
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by Polnarev [1^. As we will show, the amplitude of the polarization is de- 
termined by the energy scale of inflation, and its angular/spatial correlation 
structure is determined nearly exclusively by the expansion of the universe. 
While the energy scale (and, quite frankly, even the existence) of the inflaton 
is unknown, the post-inflation expansion history of the universe is extremely 
well understood. This is quite fortuitous for experimentalists hunting for 
the GWB as it dramatically restricts the range of our prey! 

The separation between the inflationary and standard hot Big Bang de- 
pendencies is yet another manifestation of the interplay between inflation's 
quantum mechanical aspects and the Big Bang cosmology's classical dy- 
namics. Although the inflationary perturbations are quantum mechanical 
in origin, they are small enough to be treated using linearized classical gen- 
eral relativity (the so-called WKB semi-classical approximation). So while 
the tensor-to-scalar ratio will probe quantum cosmology, the (classical) evo- 
lution of the scale factor allows for a precise prediction of the GWB's angular 
correlation imprint on CMB polarization. This separability, into a classi- 
cal part (sensitive to the background evolution of spacetime) and a small, 
perturbative quantum component, makes the CMB's curl-mode polarization 
the most robust probe of inflation. 

1.2.3 Observations and challenges 

The inflationary model has revolutionized cosmology. Inflation solves the 
"horizon problem" — reconciling observations that show that regions of the 
universe have identical CMB temperatures (to a part in 10^) by providing 
a causal mechanism for these regions to attain thermal equilibrium 380,000 
years after the Big Bang. Inflation solves the horizon problem via an expo- 
nential, accelerating expansion of the universe at early times, prior to the 
"ordinary" Hubble-Friedmann expansion observed today. This rendered the 
entire observable universe in causal contact initially, and also accounts for 
the seemingly finely tuned spatial flatness of the universe observed by CMB 
temperature anisotropy experiments [201 Ell ESI ESl EH EH] . 

Inflation also predicts a nearly scale-invariant spectrum of scalar perturba- 
tions. That any initial perturbations remain after the universe expanded by 
a factor of ~ e^^ is astonishing! Yet, surprisingly, the fluctuation level at the 
surface of last scattering arises naturally [14J in inflation as a consequence 
of parametric amplification, see Section [1.4.11 The residual fluctuations are 
observable in the CMB and indicate the epoch of inflation and the amount of 
expansion (the duration of inflation). This is inflation's solution to cosmol- 
ogy's "smoothness problem," accounting for the small, but non- vanishing. 
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level of perturbations. Regrettably, neither flatness nor smoothness are 
unique to inflation. Both have long histories, predating inflation. Flat- 
ness was anticipated on quasi-anthropic principles [26], and the universe's 
near-smoothness was predicted as the primordial matter power spectrum 
[271 l28l [29] . Recent CMB and galaxy cluster measurements [301 EI] have de- 
tected perturbations (possibly) resulting from phase-synchronized "quantum 
noise" (zero-point oscillations in the inflaton). These scalar, mass/energy 
perturbations, combined with the universe's spatial flatness, increase infla- 
tion's credibility since the e^'^-fold expansion producing flatness should have 
also destroyed all initial perturbations. However, while there is abundant 
circumstantial evidence, there is one unique prediction of inflation: the pri- 
mordial GWB, which produces an unmistakable imprint on the polarization 
of the CMB. 

1.3 CMB polarization 

The CMB is specifled by three characteristics: its spectrum, the spatial 
distribution of its intensity (or temperature anisotropy), and the spatial dis- 
tribution of its polarization. All three properties depend on fundamental 
cosmological parameters. Additionally, since CMB photons travel through 
evolving structures in the early universe on their way to our telescopes to- 
day, the CMB is also a probe of cosmic structures along the line of sight, 
which are, in some sense, "foregrounds" — either emitting, attenuating, or 
distorting the spatial and frequency power spectra of the background. 

Originally proposed by Rees [32] as a consequence of an anisotropically 
expanding universe, the polarization of the CMB was unobserved for many 
decades. Although Rees' original model was found to be untenable, it was 
later corrected by Basko and Polnarev [33] in 1980. Nevertheless Rees' explo- 
ration attracted the attention of experimentalists |341 l35l |36] who initiated 
observations to measure CMB polarization. The polarization of the CMB, 
and its correlation with temperature anisotropy, was flrst detected by DASI 
|37] . The race to discover the wispy imprint of the GWB was on! 

1.3.1 Temperature anisotropy produced by the GWB 

When electrons in the primordial plasma prior to decoupling were irradi- 
ated with CMB photons, polarization of the microwave background was 
inevitable. Thomson scattering (low energy Compton scattering) produces 
polarization whenever photons from an anisotropic radiation fleld scatter 
off unbound electrons. Anisotropy in the CMB radiation field was pro- 
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duced by either mass/energy perturbations (over- and under-dense regions) 
or gravitational waves. When the photon field is decomposed into spher- 
ical harmonics, these two types of perturbation produce anisotropy of the 
quadrupolar variety (i^^m with 1 = 2). There are five harmonics with i = 2, 
but only one of these, with m = 0, is azimuthally symmetric. The ^2,™ with 
m = ±2 indicate that gravitational waves are spin-2 objects [38]. Gravita- 
tional waves "shear" spacetime and produce local violations of refiection, or 
parity, symmetry in the CMB polarization field |19l [55] . 

The GWB produces CMB temperature anisotropy as well. However, the 
temperature anisotropy is a scalar field on the celestial sphere and is dom- 
inated by the acoustic oscillations of radiation and matter, overwhelming 
the minute temperature anisotropy produced by gravitational waves. The 
temperature anisotropy induced by the GWB is also degenerate with other 
cosmological parameters SO] and essentially undetectable at levels below 
the current WMAP3 limits, due to cosmic variance j41] . 

1.3.2 Polarization anisotropy produced by the GWB 

Fortunately, however, the CMB polarization's tensorial nature breaks the 
parameter degeneracy. Using an analog of Helmholtz's vector calculus the- 
orem valid for spin-2 fields on the celestial sphere, CMB polarization maps, 
like Figure [TTT] can be decomposed into two scalar fields or "modes" [39] HQ]. 
The advantage of manipulating two scalar fields, as opposed to one tensor 
field, is self-evident. 

One of the scalar fields is, essentially, the gradient of a scalar-potential 
and is known as "E-mode," or "gradient-mode," polarization by analogy 
to the electric field. The E-mode polarization is invariant under parity 
transformations. The second component, called "B-mode," or "curl-mode," 
polarization is analogous to the the curl of a vector-potential. If infiation 
produced a sufficient amount of gravitational radiation, then future maps 
of CMB polarization will be admixtures of both modes (though the E-mode 
polarization will dominate by at least a factor of ten). For reference, simple 
one-dimensional maps of pure E- and B-modes are shown in Figure [L2l Ap- 
praising the behavior of the circular maps with respect to refiections across 
the map's diameter reveals the symmetry of the underlying polarization 
mode. 

The pioneering work by Polnarev [19] was the first to identify a unique 
observational signature of gravitational waves; one that would only be man- 
ifest in the polarization of the CMB. Polnarev's key insight was to recognize 
that asymmetric shear induced by gravitational waves would induce a polar- 
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Fig. 1.1. A simulated noiseless map of CMB polarization and temperature 
anisotropy. The simulation represents an 18° x 18° map of the CMB's tempera- 
ture, E-mode, and B-mode polarization. The temperature (gray scale) and E-mode 
polarization reveal the classical cosmological parameters such as the mass density, 
geometric curvature, and composition of the universe (i.e. "dark" versus ordinary 
matter). The B-mode, or "curl," polarization is only generated by primordial grav- 
itational waves and is indicated by regions where reflection symmetry is locally 
violated. The vertical scale bar at the lower left indicates polarization at the 4 /iK 
level. The polarization vectors are the sum of E- and B-mode polarization, but 
are dominated by E-mode polarization (B-mode polarization is less than 0.1 /LtK in 
this simulation, corresponding to a tensor-to-scalar ratio r 0.1). Figure credit: 
Nathan Miller. 



ization pattern significantly different from that produced by scalar pertur- 
bations, such as those shown in the right-hand side of Figure [L2l Polnarev 
predicted that gravitational waves would be the only plausible source of 
parity violation on a cosmological scale. In the more modern language of 
E- and B-modes, this is equivalent to predicting the existence of B-modes 
imprinted on the CMB by a primordial GWB. 
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Fig. 1.2. Parity symmetric E- mode (or "gradient-mode") polarization patterns (left 
side), and B-modc (or "curl- mode") patterns (right side) in real space. Figure 
credit: Nathan Miller. 



1.4 The origin of the gravitational wave background 

How can any perturbations originating from quantum fluctuations in the 
primordial inflaton field survive the explosive expansion by a factor of 
After all, a hallmark of inflation is that this expansion dilutes the curvature 
of the universe from any primordial value to precisely flat. In fact, since 
the GWB is a radiation background exactly like the CMB, the subsequent 
expansion following inflation dilutes the GWB's energy density by a factor of 
a^, where a is the cosmic scale factor. Since the GWB energy density today 
is at least one billion times smaller than the CMB's energy density, this 
means that it was utterly insignificant for all times, including at decoupling. 

This brings us to another potentially troubling "why" question: Why 
are gravitational waves expected to persist to last-scattering, and leave a 
detectable imprint on the CMB, if the very cosmological model predicting 
them produces an expansion that should render them negligible? A hint 
at the answer to this question comes from a rather unlikely source: a play- 
ground swingset! 
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1.4 -1 Parametric resonance and amplification 

Parametric amplification is exemplified by an undamped pendulum of length 
L whose suspension point, y, is driven vertically y = A cos Inft. To highlight 
the connection between (1) the GWB, (2) the vertically driven pendulum, 
and (3) the laser, we term the periodic driving force the "pump." Pump 
energy at frequency / drives the pendulum into resonance and can even am- 
plify small random thermal vibrations of an (initially stationary) pendulum 
into oscillation. 

For simplicity, the angle between the pendulum and the vertical, <j), is 
taken to be small, and the equation of motion for (j) becomes: 

4>+ {ujl + n'^COs2TTft)<j) = (1.1) 

where = {2T:ffA/L. 

Equation 11.11 is known as Mathieu's equation. In Mathieu's equation, ujq 
is the natural frequency of the pendulum, and the parameter $7 determines 
the resonant behavior of the pendulum, leading to parametric resonance. 
While the equations of motion are linear with respect to (p, the effect of 
the pump is not additive (as it would be if the suspension point were to be 
horizontally modulated), but rather multiplicative. 

To find periodic solutions to Equation 11.11 we construct the following 
ansatz: 

4> = ^+e+'^f^ + ^.e-'^f^ (1.2) 

In general, both stable and unstable solutions of Equation 11.11 can be 
obtained. So called "resonance bands" are separated by regions of stability 
where the amplitude of (f) is constant. Unstable solutions exponentially 
diverge (as ~ e*). Both types of solutions will be important in the context 
of gravitational waves. Solving Equation 11.11 using Equation 11.21 leads to 

- (vrV')]^± + %Vt = ^ (1-3) 

where third-harmonic generation effects have been ignored. Stable solu- 
tions to the (two) Equations 11.31 are non-trivial and solvable for ip± when 
the following self-consistency relation holds between the pump amplitude, 
frequency, and natural frequency of the pendulum: 

f]2 = 2[a;2_vr2/2] (1.4) 
The first, or fundamental, resonance condition for the pump frequency 
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is / = loq/tt, which imphes that pumping at twice the pendulum's natural 
frequency defines the boundary of incipient instability even if the pump 
amplitude is small 7^ even as rt,A^ 0). 

The pumping strategy mentioned above shares similar features with a 
similar application on the playground. Assisted swinging on a swingset is a 
resonant system with pump-power supplied by two assistants, one at each of 
the two displacement maxima; that is, with pump frequency f = 2^. When 
this condition holds, J7 = and only small amounts of pump energy are 
required to maintain the swing's oscillatory behavior, even in the presence 
of significant frictional damping (which has been ignored here). 

In fact, while less social (and more dangerous) than assisted pumping 
(using two friends), vertical pumping employing parametric resonance al- 
lows the rider to initiate resonance by themselves (by raising and lower- 
ing their center of mass — alternately standing and squatting on the swing) . 
Surprisingly, the amplification of small initial perturbations via parametric 
resonance provides a fruitful analogy for the theory of gravitational wave 
amplification. 

1.4-2 Parametric amplification of the GWB 

Gravitational waves have unique and fascinating cosmological properties. 
While the contribution of these waves to the energy density of the universe 
today is minuscule, parametric amplification of these primordial quantum 
fiuctuations of the inflaton field causes the waves to grow large enough to 
become potentially observable. If the imprint of these primordial perturba- 
tions is observable, an understanding of the parametric amplification pro- 
cess allows us to optimize our observational requirements — regardless of the 
magnitude of the inflationary GWB. 

In our simplified cosmology, spacetime is smooth and flat. On top of this 
background a tensor field, representing the GWB, is suffused. The metric 
of this spacetime is obtained by solving the Einstein equations 

Gap = 8TrGTaf3 

where G is Newton's constant of universal gravitation. For empty space the 
stress-energy tensor Tap = 0, leading to 

ds^ = a{r]f [dTf - {5c,p + hap)dx''dx^] (1.5) 

where 5ap is the Kronecker delta function, and ry, the conformal time, is 
related to the (time-dependent) cosmological scale factor a via drj = dt/a. 
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The (linearized) perturbation tensor is both transverse-symmetric {ha/3 = 
hf^a) and traceless {J^ahaa = 0). We seek solutions, which are separable 
into a tensorial part and a scalar part, of the following form 



ha/3 = V8^-ea/3e"''^ (1.6) 

where e^/j is the gravitational wave polarization tensor. The rank-two ten- 
sors hai3 and e^/? are transverse and traceless, leading to two independent po- 
larization modes denoted "+" and "x." Solving Einstein's equations yields 
wave equations for the (scalar) amplitude of the set of equations hap 



u + {k^ - 'd/a)u = 
i) + {k^ -Ueff)l^ = (1.7) 

Here, overdots denote derivatives with respect to conformal time (e.g. 
X = a^), and in the second equation we have replaced U^ff = f , which 
acts as a time-dependent effective potential [12]. We recognize Equation 
11.71 as a version of Mathieu's equation for the vertically driven pendulum, 
Equation ll.il once the following substitution is made in Equation ll.il 



—Q^ cos2tt ft = Ueff{a) 

In contrast to the vertically driven pendulum, variation of the pump parame- 
ter Uf>ff does not lead to runaway growth. Rather, the time- varying effective 
potential amplifies long-wavelength oscillations relative to short wavelength 
oscillations. 

For an isotropic, homogenous universe consisting of a fluid with pressure 
p and density p we can express 



Ueffia) 



da 



d 

It V dt 
da^^ 

2 



+ a 
SttGp 



(fa _ 

dt"^ ~ 
AttG 



a^H^ + a 



Sa 
dt^ 



(p + 3p/c2 



(1.^ 



Here, we have employed the definition of the Hubble parameter H(a) = 
= y^87rG/f3(a)/3. For convenience, the equation-of-state relating pres- 
sure, p, and density, p, is taken as p = 7pc^ (where 7 is a scalar that depends 
on the cosmological epoch under consideration). As the universe expands. 
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it dilutes: p = P*i-^)^^^^^'^^ ■ This equation is valid at any epoch, or corre- 
spondingly, for any value of a. When we consider a specific epoch we label 
it a* . From Equation 11.81 we obtain 



eff\ 



1 + 37)1 4^Gp*(l -37)a2 



2-3(1+7) 



which can be written as 



4^Gp.(l-37)a^ / a ^ -{1+37) 
Ueff{a) = [-) (1.9) 

The evolution of the effective potential depends crucially on the cosmo- 
logical epoch, via the relationship between density and pressure. For ex- 
ample, during radiation domination, 7 = 1/3 and the effective potential 
Ueff (a) = 0. Recalling Equation 1 1.71 when either Ueff{a) = or fc^ » ^e//; 
the solutions for v, the gravitational wave amplitude, are simple plane- waves. 
These purely oscillatory solutions prevail whenever k 3> Uef / , but especially 
during radiation domination when the dilution of the GWB is identical to 
that of any radiation background, such as the CMB. 

For future reference, using the definition of the Hubble parameter H{a) 
in terms of density p{a), we can write 

where = a^:H{a^:). 

Using Equations 11.61 and 11.91 we can solve Equation 11.71 for the gravita- 
tional wave amplitude 

/i = VS^— e[^'=(''-''*)] (1.11) 
a 

Gravitational waves "enter" the horizon when k ~ 1/?], and as the universe 
Hubble-expands, they are damped by the adiabatic factor 1/a, just like 
radiation. Example solutions of Equation 11.111 are shown in Figure 11.31 



1.4-3 The effective potential 

We have shown that the equation-of-state parameter 7 determines the effec- 
tive potential, subsequently determining the evolution of the GWB, at least 
during radiation domination. In this subsection we examine solutions in the 
other important cosmological epochs. We will see that specifying 7 versus 
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Fig. 1.3. Gravitational wave amplitude as a function of conformal time for four 
different wavenumbers. Gravitational waves are constant, and equal in amplitude 
independent of wavelength, before entering the horizon when fcry ~ 1, leading to de- 
cay. Equivalently, waves with k <C Ueff{a) (long wavelength waves) experience the 
effective potential, forestalling their decay. Short wavelength waves never interact 
with the potential, continuously decaying adiabatically instead. Long wavelength 
modes essentially "tunnel" through the barrier with no diminution of their ampli- 
tude until after radiation-matter equality. Figure credit: Nathan Miller. 

cosmological scale not only allows us to predict the advantages of indirect 
detection of the GWB (using B-modes) over direct detection methods, but 
also to optimize experimental CMB polarization surveys themselves. 

In the following, the subscripts "-" and will denote quantities before 
and after inflation ends, respectively. If aend denotes the scale factor at the 
end of inflation, this means that for a < Oend,! < — |- After inflation ends 

> O'end and 7+ > — |. For example, 7+ = 1/3 corresponds to the equation- 
of-state for radiation; that is, when inflation ends the universe is radiation 
dominated. 

More generally we can say that at the moment when a = a^ndi accelerated 
expansion [dP'a/dt^ > 0) changes to decelerating expansion (cPa/dt'^ < 0), 
and during inflation the cosmological horizon decreases [Hj. Hence, from 
Equation 11.101 for a < aend, 




(1.12) 



which increases with a, and when a > a^nd 
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(l-37+)A:,^ /Q\-|i+37+l 
2 \aj 



(1.13) 



which decreases with a. Since the transition from 7_ to 7+ occurs quasi- 
instantaneously, the effective potential is discontinuous; 7+ > 7-- For more 
accurate results we should properly treat the reheating phase, when the 
inflaton is converted into particles and radiation. Interestingly, parametric 
resonance techniques can also be used to describe reheating 03]. 

Now it should be clear why | was called an effective potential. Since 
Ueff is maximized when inflation ends, there are two epochs, a_ < a^nd and 
a_|_ = a^ndj where k'^ = U(.ff < Umax- Thus there are two wavelength regimes 
that determine the form of solutions to Equation 11.71 High-frequency waves 
with k/kn ^ r]eq enter the horizon well before matter-radiation equality, 
then decay as the universe expands as in Equation 11.111 \M[ B5] . 

Long wavelength waves, on the other hand, satisfy 



To solve Equation II . 141 we must first determine the behavior of the effective 
potential in a form that is valid during any cosmological epoch. 



To analyze the impact of the effective potential on long wavelength gravi- 
tational waves we must solve for Ueff ~ a/a, recalling that the derivatives 
are with respect to conformal time, r]. We therefore require the relationship 
between the scale factor and conformal time. For reference, we recall that 
during radiation domination Ueff = for all wavelengths. 



During the epoch of matter domination, it is convenient to parameterize the 
evolution of the scale factor versus time as a ~ t"^ with a = 2/3. Using the 

a. 

definition of conformal time, we have drj = dt/a, implying that a = ryi-^ or 
a{r]) = rf during matter domination. Using this, we find that the effective 
potential during matter domination decays as Ueff {a) oc 1/a, as described 
quantitatively in Equation 11.131 and displayed in Figure 11.41 

Finally we must calculate the effective potential during inflation, when 
the scale factor grows as a = e^*, implying dr] = e~^^dt. Integrating, we 
find ?? = (1 - e~^^)/H or that 




(1.14) 



1.4 -4 Timing is everything: the cosmic chronology 
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1 



(1.15) 



a ~ 
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Fig. 1.4. The parametric amplification of primordial gravitational waves is governed 
by the evolution of the effective potential as a function of cosmic scale factor, a. 
This figure, adapted from [42], shows the effective potential in three important 
cosmological epochs. Here the wave's frequency {vq = ck) is expressed in present- 
day units, when a — 1. Long wavelength waves (fc^ < Ueff) "tunnel" through the 
potential and remain constant until inflation ends. Short wavelength waves never 
experience the potential and instead decay adiabatically (as 1 /a) . During radiation 
domination, the effective potential vanishes and all waves inside the horizon decay. 
Finally, any waves that survive until matter-radiation equality imprint the CMB 
sky prior decoupling. Therefore, these waves are comparable to, or larger than, the 
horizon at decoupling, subtending an angle of ~ 2° on the sky today. Figure credit: 
Nathan Miller. 

during inflation, leading to 

Ueff{a) = 2H^a^ (1.16) 
That is, the effective potential grows quadratically. 

1.44.2 Evolution of the GWB 

With the effective potential expressed in terms of conformal time we can 
easily solve the gravitational wave equation. During inflation, a ~ 1/ry 
so Uf,f f (X (as anticipated from Equation I1.12p . For long wavelength 

waves, the parametric equation is 
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i}-^u = (1.17) 

Solutions to Equation 11.171 are easily verified to be ~ 1/??, or i'{a) ~ 
a. Since h{a) = VSirGu/a (Equation II. 6p we see that the gravitational 
wave amplitude during the inflationary epoch is constant for waves with 
wavelengths smaller than Ueff{a < aend), or equivalently 

h(k,rj) = ^^ (1.18) 
rja 

During radiation domination, U^ff = and the gravitational waves are 
oscillatory functions that redshift and dilute adiabatically as 

a 

During matter domination, although the effective potential has the same 
dependence on r/ as it does during inflation, the scale factor depends on 
conformal time in a different way. This leads to a different gravitational 
wave solution. During matter domination, r] = yfa and so h = u/a ^ 
1 /a^/^ , implying that gravitational waves decay even when their wavelengths 
are long {k < Ueff{aeq)). Once again, however, longer wavelengths are 
preferentially preserved relative to short wavelength modes. 

Thus we have the following cosmic chronology. During inflation, long 
wavelength gravitational waves are "frozen" outside the horizon, or equiva- 
lently, are spared from adiabatic decay as they tunnel through the effective 
potential, which grows quadratically as the universe inflates. At the begin- 
ning of radiation domination (end of inflation; Ognd), the effective potential 
drops to zero. All waves that are within the horizon decay adiabatically. 
The largest waves persist until decoupling, when the universe is matter dom- 
inated and the potential again becomes critically important. The primordial 
scale-invariant distribution of waves is transformed twice; first during infla- 
tion, and later during matter domination. In both cases, longer wavelength 
perturbations are overpopulated with respect to their short wavelength coun- 
terparts. 



1.4 -5 Quantum gravitational wave effects 

We can make an analogy between the GWB and the laser here by rewriting 
the GWB parametric equation (Equation II. 7p as i> + {k'^ — U(.ff{a))u = 
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i> + ^eff^ ~ ^- ^eff ^ *-* solutions are constant oscillatory functions, 
while for /c^jj < we have decaying solutions. Similar derivations [441 
145) . using creation and annihilation operators, are particularly useful for 
exploring the connections between the quantum properties of the GWB 
and the laser (which can also be derived using creation and annihilation 
operators). This derivation is possible because the gravitational waves, or 
gravitons, are bosons, as are the laser's photons. In both cases, the pump 
need not be periodic. 

However, the analogy cannot be taken much further. For the prototypical 
three-level laser, coherent amplification is obtained via stimulated emission 
from a metastable state. The metastable state's occupation number is in- 
verted with respect to the ground state. For the GWB there is no such 
population inversion, nor stimulated emission, and thus the GWB is an 
incoherent, stochastic radiation background, not unlike the CMB. 

We are interested not only in the properties of a single gravitational wave, 
but also in the behavior of a stochastic ensemble of waves, tracing its origin 
to "quantum noise" in the inflaton. To quantify the preferential population 
of long wavelength gravitational waves, we calculate the initial number of 
gravitons, or occupation number, for given conformal wave number k when 
a <^ ttend, denoted as Nin- We construct the correlation function of this 
stochastic background by considering the energy of a collection of os- 
cillators, each with energy E = huj{n + ^) in state n. Using terminology 
familiar from laser physics, the number of gravitons [36] in a volume V{a) is 
the (renormalized) energy divided by the product of the frequency and the 
reduced Planck constant, h: 



~ 52G7Thuj{a)\y dt jy dt )/ ^^■^^> 

In Equation 11.191 * indicates complex conjugation, the angled brackets 
(. . .) denote averages over polarization, propagation angles, time, and vol- 
ume. The average results in a constant factor that is subsequently incorpo- 
rated along with the physical constants into a constant k. Finally, we define 
Lo = k/a, thus simplifying the time derivatives: = cohp and V{a) = Voa^. 

As result we have 



h-.a. 



2/A:\2//cN-i „ 
- - Voa' 
a/ \a/ 



nh^cLLkVf) 



(1.20) 
fl.21) 
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which does not depend on a. Similarly, the final mode occupation number 
at the end of inflation, Nout, when a ^> Cend ^-nd k'^ S> U^jf is 

Nout ~ nh\a\kVo 

which is also independent of time and a. For a_ < a < a+, when k <C 
Uef f (long wavelength gravitational waves) , the solution for the gravitational 
wave amplitude is ~ a. For a scale invariant spectrum h is constant, 
meaning that gravitons are created with the same amplitude ~ ^+)- 
This leads to the amplification factor: 

Nin ^ Va_/ 

which is larger than unity if fc^ < Umax- Therefore, long wavelength grav- 
itational waves are amplified with respect to their short wavelength coun- 
terparts (which experience the potential for a much shorter time). This is 
evident from Figure 11.41 Intermediate wavelength waves [k"^ ^ Umax) do 
not experience the effective potential until much later than their long wave- 
length counterparts, and very short wavelength waves, with k'^ > Umax, do 
not enter the potential at all. 

As such, if all waves start with the same initial amplitude (Harrison- 
Zeldovich spectrum), short wavelength waves will decay by a quadratic fac- 
tor relative to longer waves. We note that the change in scale factor during 
inflation can be enormous; expansion by factors of 10'^'^ are common for mod- 
els that produce sufficient inflation to make spacetime flat. In the context of 
the GWB, only the difference in the expansion of the universe between the 
horizon-entry time for short wavelength modes compared to that for long 
wavelength modes is relevant. 



(1.22) 



1.5 Observational consequences 

To illustrate the significance of parametric resonance, it is useful to consider 
the following generic inflation scenario where the equation-of-state param- 
eter changes from 7_ = —1 to 7+ = 1/3; that is, radiation domination 
immediately follows inflation. Using Equation 11.101 iii such a scenario a 
wave with k^nd — cbendHend enters the effective potential just as inflation 
ends. As such it is not amplified. A long wavelength wave enters the poten- 
tial earlier when a_ ~ H/k. We are free to define the scale factor at the end 
of inflation as Ognd = 1- Thus from Equation 11.221 with = ag^^ = 1, the 
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amplification scales with wavevector as ^ ~ ^-fr^, assuming H is constant 
during inflation. The implication of this equation is clear: if the primordial 
gravitational wave spectrum is scale invariant ~ ^+), then post-inflation 
it is transformed into a strongly wavelength-dependent spectrum. 

Not all of these amplified waves survive to imprint the CMB with B-mode 
polarization. Waves that are larger than the present-day horizon are frozen 
and are not amplified. Additionally, all sub-horizon waves decay by the 
adiabatic factor from the onset of radiation domination at reheating to 
the epoch of matter-radiation equality. 



1.5.1 Tensor power spectrum: amplitude and angular structure 

More generally, we are interested in the variance of h[k, rf) as a function of 
wave number, also called the tensor power spectrum, defined as 



m) 



From Equation II . 181 we recall that h{k,rj) = . Since r] ~ 1/aH (Equa- 
tion [TTT5]), we find 

Pt{k) = — p — (1.23) 

Thus, measuring the tensor power spectrum probes the Hubble constant 
during inflation and is proportional to the energy density of the inflaton 
(since oc p). The amplitude of the tensor power spectrum is directly 
revealed by the amplitude of the CMB curl-mode polarization power spec- 
trum. If inflation occurs at energy scales comparable to Egut, CMB B-mode 
polarization allows us to test physics at scales one trillion times higher than 
the highest energy produced in terrestrial particle accelerators! 

We have shown how the tensor power spectrum's amplitude depends on 
the energy scale of inflation. Now our aim is to qualitatively determine 
the angular correlation properties, or shape, of the tensor power spectrum 
Pt{k). Doing so allows us to optimize observations of the GWB, for as we 
will see, Pt{k) has a well-defined maximum as a function of wavenumber, or 
equivalently, angle subtended on the sky. 

In the context of the driven undamped pendulum, the ansatz solutions. 
Equation II. 2^ can grow exponentially or remain constant. The long wave- 
length gravitational waves amplitudes are also stable when they are larger 
than the cosmological horizon. Short wavelength gravitational waves decay 
adiabatically because they are always within the horizon during inflation. In 
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the context of parametric amplification, long wavelength waves (A;^ < C^e//) 
are said to be "super-horizon" or non-adiabatic. 

Since the extremely high-frequency waves are always above the potential 
barrier, they continuously decay and leave no observable signature. The 
longest waves enter the horizon near decoupling when the effective potential 
is decaying. Waves that experience the matter-dominated effective potential 
close to decoupling are amplified the most, leading to a well-defined peak 
in the tensor power spectrum on scales comparable to the horizon at decou- 
pling. This scale subtends an angle of ~ 1° to 2° on the CMB sky today, and 
therefore this is the characteristic scale of the tensor and B-mode angular 
correlation function. 



1.5.2 Direct versus indirect detection of the primordial G WB 

It is perhaps instructive to ask whether the primordial GWB could be di- 
rectly detected today. While scalar, or mass-energy, perturbations are am- 
plified by gravitational condensation, the tensor GWB is not. Directly de- 
tecting the GWB today (redshift z = 0) by, for example, LIGO, would be 
extremely difficult since, like the primordial photon background (the CMB), 
the energy density of the primordial GWB dilutes (redshifts) by the fourth 
power of the scale factor as the universe expands. However, the GWB im- 
prints curl-mode polarization on the CMB at the surface of last scattering 
(at z ~ 1,100 or 380,000 years after the Big Bang). Therefore, the energy 
density of the GWB at last scattering was more than one trillion times 
(1, lOO'') larger than it is now. 

Ultimately, direct detection experiments such as ESA's LISA, NASA's Big 
Bang Observer, or Japan's DECIGO experiment will provide further tests 
of the infiationary model, such as measuring the GWB power spectrum at 
wavelengths approximately twenty orders of magnitude smaller than those 
probed by CMB polarization [Ul SSI SS]. As we have seen, the current 
spectral density of these short wavelength waves is extremely small, and 
these experiments are fraught with contamination from "local" (i.e. non- 
cosmological) sources. However, given a detection of the primordial GWB at 
the surface of last scattering using CMB polarization, these direct detection 
campaigns will measure the fine details of the infiaton potential, making 
them at least well justified, if not mandated. 
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1.5.3 Indirect detection of the GWB: optimizing CMB 
polarization observations 

Parametric amplification had important ramifications for the design of the 
first experiment dedicated to measuring the GWB — BICEP (Background 
Imaging of Cosmic Extragalactic Polarization). While the amplitude of the 
GWB is unknown, parametric amplification allows the structure of the B- 
mode's angular correlation function to be accurately calculated given only 
modest assumptions about the cosmic equation-of-state and evolution of 
the scale factor. This allowed the BICEP team to optimize the angular 
resolution of our search for the B-mode signature; motivating both BICEP's 
optical design as well as its survey design (required sky coverage). 

BICEP probes the inflationary GWB primarily at large angular scales 
corresponding to the largest wavelength waves that entered the horizon near 
decoupling (~ 2°). Since the tensor angular correlation function peaks on 
these scales, to obtain statistical confidence in our measurement it is only 
necessary for BICEP to probe a small fraction of the sky (~ 3%), rather than 
diluting our observing time over the entire sky. This allows us to target the 
cleanest regions of the sky — those with minimal contamination from galactic 
dust or synchrotron radiation, both of which are known to be polarized. 

1.5.4 Detectability of the CMB B-mode polarization 

Nearly thirty years passed between the discovery [50] of the CMB by Robert 
Wilson and Arno Penzias (Charles Townes' Ph.D. student) and the first 
detection of temperature anisotropy [JT] by COBE at the ten parts-per- 
million level. If inflation occurred at the GUT-scale, it would produce curl- 
mode polarization at the ten parts-per-6i//ion level. 

Thanks to the innovative technologies produced by our collaboration, de- 
tection of this signal is conceivable. BICEP — the ground-based polarimeter 
we built — will achieve higher sensitivity to the GWB than either NASA's 
WMAP (Wilkinson Microwave Anisotropy Probe) or the Planck satellite. 

BICEP is both a pioneering experiment and a long-range, two-phase cam- 
paign designed to mine the CMB sky using innovative technology. BICEP 
is an attempt to probe even farther back than the last scattering surface: 
to the very beginning of the universe, the inflationary epoch. Detecting 
the GWB requires ultrasensitive technology, only recently invented. An un- 
derstanding of parametric amplification allows us to precisely estimate the 
spatial power spectrum imprinted on the CMB's polarization by the pri- 
mordial GWB. This leads to a very general optimization of experimental 
campaigns [52]. To detect the GWB's polarization imprint, only modest an- 
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gular resolution is required — corresponding to a small refractor. This small 
refractor can probe the inflationary B-mode polarization nearly as well as a 
reflecting telescope twenty-times larger in diameter! 

As shown in the following section, BICEP's small size has several impor- 
tant ancillary benefits — most notably that its smaller aperture results in a 
much higher-fidelity optical system, one with no obscuration or secondary 
mirror to induce spurious polarization. Additionally, the refractor is easy 
to shield from stray light — which is particularly important when probing 
signals one billion times smaller than the background. 

1.6 Experimental quantum cosmology: the BICEP project 

BICEP [531 [51] is a bold first step toward revealing the GWB. Ultimately, 
only a small telescope like BICEP can be cooled entirely to nearly the tem- 
perature of the CMB itself — a condition not previously achieved, even in 
space. BICEP's elegant design (see Figure [T3|) has proven extremely attrac- 
tive for proposed future experiments [55|156J, which have receiver concepts 
closely resembling BICEP. 

BICEP is the first experiment to directly probe for the primordial GWB. 
Even Planck, when it is launched in 2008, will not be as sensitive to the 
GWB signal as BICEP (which will already have completed first-phase ob- 
servations). A comparison of BICEP's capability to detect the inflation- 
ary GWB with that of WMAP and Planck is demonstrated in Hivon & 
Kamionkowski [57] . They show that BICEP will achieve higher sensitivity 
to the GWB than these spaceborne experiments. BICEP's sensitivity re- 
sults from advances in detector technology and from the ability to target 
only the cleanest regions of the microwave sky, rather than spreading out 
limited integration time over the full sky. 

1.6.1 Optics 

As outlined above, only modest angular resolution is required to detect the 
GWB's polarization signature. BICEP was designed to map ~ 3% of the sky 
with 0.9° resolution (at 100 GHz), 0.7° resolution (at 150 GHz), and 0.5° 
resolution (at 220 GHz). Unlike Planck or WMAP, BICEP was designed 
specifically for CMB polarimetry, able to modulate the polarization signal 
independent from the temperature signal with high fidelity. BICEP achieves 
this fidelity by virtue of an elegant optical design: a 4 kelvin refractor (Fig- 
ure ll.5p . Millimeter- wave radiation enters the instrument through a 30 cm 
diameter vacuum window and passes through heat-blocking filters cooled 
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Fig. 1.5. A cross-sectional view of the BICEP receiver, which comprises a refracting 
telescope and 49 polarization-sensitive bolometers. The optics and focal plane are 
housed within a cryostat, which is placed on a three-axis mount. Figure credit: 
Thomas Renbarger. 



to 4 K by liquid helium. Cold refractive optics produce diffraction-limited 
resolution over the entire 18° field-of-view. 
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1.6.2 Detector system 

Each one of BICEP's 49 pixels comprises a complete polarimeter (optics, po- 
larization modulator, analyzer, and detector). Each pixel uses three corru- 
gated feedhorns feeding a polarization-sensitive bolometer pair (PSB), which 
simultaneously analyzes and detects linearly polarized light. The FSB's in- 
genious design is further described in [58]. In addition to their use in BICEP, 
PSBs have been successfully used in the BOOMERANG [59] and QUAD ex- 
periments to measure the E-mode polarization signal, and they will be used 
on Planck. 

BICEP's PSBs use two absorbing grids, each coupled to a single linear 
polarization state. The (temperature-dependent) resistance of a semicon- 
ducting thermistor (neutron transmutation doped germanium; NTD-Ge), 
located at the edge of the absorber, detects CMB photons via a resistance 
change. BICEP's PSBs have astounding sensitivity: in one second, each 
PSB can detect temperature fluctuations as small as ^ 450 //K. 

1.6.3 Polarization modulation 

The faintness of the GWB polarization signal demands exquisite control of 
instrumental offsets. There are two ways to mitigate offsets: (1) minimize 
the offset and (2) modulate the signal before detection (Dicke switch) faster 
than the offset fluctuates. BICEP does both. A bridge-circuit differences 
the two PSBs within a single feed, producing a (first) difference signal that 
is null for an unpolarized input. This minimizes the offset. For six of the 
49 spatial pixels, the polarized signal input is rapidly modulated (second 
difference) by Faraday Rotation Modulators (FRM) that rotate the plane 
of linear polarization of the incoming radiation. The FRMs make use of 
the Faraday Effect in a magnetized dielectric. Polarization modulation al- 
lows the polarized component of the CMB to be varied independently of the 
temperature signal, allowing the response of the telescope to remain fixed 
with respect to the (cold) sky and (warm) ground. This two-level differ- 
encing scheme allows for two levels of phase-sensitive detection, allowing 
optical systematic effects, associated with the telescope's antenna response 
pattern (leaking the much-larger CMB temperature signal to spurious CMB 
polarization), to be distinguished from true CMB polarization. 

The FRMs represent a significant advance in the technology of CMB po- 
larization modulation. Early CMB polarimeters (including Penzias and Wil- 
son's, which was polarization-sensitive) used rotation of the entire telescope 
to modulate CMB polarization. These experiments [351 l36l [60] rotated hun- 
dreds or thousands of kilograms, were susceptible to vibration induced mi- 
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crophonic noise, and were limited mechanically to modulation rates < 0.1 
Hz. 

The next polarization modulation innovation was a birefringent half-waveplate: 
a single crystal of anisotropic dielectric (typically quartz or sapphire) that 
phase-delays one of the two linear polarizations [341 \6T\ 162) . While the 
fragile ~ 1 kg, cryogenically-cooled crystal can be rotated at ~ 1 Hz with 
lower-vibration than rotation of the entire telescope, such a mechanism is 
prone to failure since bearing operation is a severe challenge at cryogenic 
temperatures. And since bolometers are sensitive to power dissipation at the 
lO^^^W level, even minute mechanical vibrations produced by the bearings 
are intolerable. 

Faraday Rotation Modulator^, shown in Figure 11.61 require only "rotat- 
ing" electrons (the generation of a solenoidal magnetic field in a magnetized 
dielectric) to effect polarization rotation. Therefore, FRMs reduce the rotat- 
ing mass that provides modulation by 30 orders of magnitude! Furthermore, 
these devices are capable of rotating polarized millimeter wave radiation at 
rates up to 10 kHz — faster than any conceivable time-varying temperature- 
or electronic-gain fluctuation. A superconducting NbTi solenoid wound 
around the waveguide provides the magnetic field that drives the ferrite 
into saturation, alternately parallel and anti-parallel to the propagation di- 
rection of the incoming radiation. The FRM rotates the CMB polarization 
vectors by ±45° at 1 Hz, well above 1/f-fluctuation timescales (caused by, 
for example, temperature variations). The bolometer signals from the PSBs 
are detected using lock-in amplification. 



1.6.4 Observations of galactic polarization using Faraday 

Rotation Modulators 

Initial observations of the galactic plane were obtained during the austral 
winter of 2006. Several hundred hours of data were taken with the FRMs 
biased with a 1 Hz square wave. This modulation waveform effected ±45° of 
polarization angle rotation. Other modulation waveforms can provide more 
or less rotation, as desired. 

To validate the FRM technology we targeted several bright regions of the 
galactic plane. Results from some of these observations are displayed in Fig- 
ure [L71 where we also show the same region as imaged by WMAP [18]. The 
agreement is impressive. Because BICEP's bolometers simultaneously mea- 

t us Patents Pending: "Wide Bandwidth Polarization Modulator, Switch, and Variable Atten- 
uator," US Patent and Trademark Office, Serial Number: 60/689,740 (2005). 
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Fig. 1.6. A cross-sectional view of a Faraday Rotation Modulator (FRM) pixel. 
Polarized light enters from the right, is rotated by ±45° and then analyzed (de- 
composed into orthogonal polarization components, which are detected individually 
by the polarization-sensitive bolometers, or PSBs). For each of BICEP's pixels, the 
PSBs are contained within a corrugated feedhorn, cooled to 0.25 K, located ap- 
proximately 20 cm from the FRM, which is placed in a corrugated waveguide at 
the interface between two corrugated feedhorns, placed back-to-back and cooled 
to 4 K. The (schematic) location of the PSBs in this figure serves to illustrate the 
coordinate system used as the polarization basis. Figure credit: Thomas Renbarger. 



sure polarization and temperature anisotropy, we use WMAP's temperature 
maps as a calibration source for BICEP. 
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Fig. 1.7. Map of a portion of the galactic plane made using one of BICEP's six 
Faraday Rotation Modulator (FRM) pixels (top) operating at 100 GHz, compared 
to WMAP's observations of the same region (bottom). Microwave radiation is 
polarized by dust grains preferentially aligned by the galaxy's magnetic field. The 
FRM pixels modulate only the polarized component of the emission and can fully 
characterize linear polarization without rotation of the telescope. The short lines 
indicate the magnitude and orientation of the plane of polarization, and the gray 
scale indicates the temperature scale. For comparison, a scale bar representing 
25 fiK linear polarization is shown. The galactic plane extends approximately from 
the lower left to the upper right of each map. Both maps show significant linear 
polarization orthogonal to the galactic plane, which is expected as the galaxy's 
magnetic field is oriented parallel to the plane |l8j . Similar maps were produced 
for BICEP's 150 GHz FRM pixels. We note that the BICEP FRM data were 
acquired over the course of a week of observations, whereas the WMAP data was 
obtained over three years. Figure credit: Evan Bierman. 
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1.6.5 The South Polar Observatory 

Essential to achieving maximum sensitivity to CMB B-modes is the long 
integration time afforded by the South Pole site — arguably the premiere, 
long-duration, low-background (both natural and manmade) Earth-based 
site. To exploit this location, BICEP is highly efficient, having been de- 
signed for robustness and quasi-autonomous operation, while consuming a 
minimum of liquid cryogens (a precious commodity at the South Pole). BI- 
CEP 's toroidal cryogen tanks house the instrument in a thermally uniform 
4 K environment. The 49 polarimeter pixels, optics, and sub-Kelvin re- 
frigerator are removable for easy instrument servicing. In December 2005, 
BICEP was installed at the US Amundsen-Scott South Pole Station's Dark 
Sector Laboratory, operated by the National Science Foundation. The Dark 
Sector Laboratory will house BICEP for three austral winter observing sea- 
sons and then house a future upgraded version, called BICEP-II, maximally 
leveraging the polar infrastructure investment for years to come. 

More detail about the design of BICEP can be found in Keating et al. 
[53]. Preliminary data, maps, and additional technical information from 
BICEP's first observing season can be found in Yoon et al. [54j. 

1.7 Future probes of the past: BICEP-II, an advanced CMB 

polarimeter array 

Following BICEP's initial phase, we will deploy an advanced high-density 
array to our South Pole observatory, eventually yielding ten-times better 
sensitivity than the first phase of BICEP. The GUT-scale is already nearly 
within the reach of current experiments [63j, such as BICEP, and nearly all 
detectable inflationary models will ultimately be testable with BICEP-II's 
technology. 

BICEP's NTD-Ge semiconducting detectors are background-limited, mean- 
ing the sensor's intrinsic noise is sub-dominant compared to the photon noise 
from atmospheric emission. In this regime, increasing the signal-to-noise ra- 
tio of the experiment can be accomplished only by adding detectors. After 
BICEP's three-year campaign concludes, our team will increase the number 
of detectors in BICEP's focal plane by a factor of more than five. This will 
produce a CCD-array-like focal plane which we call "BICEP-II." 

BICEP-II will use the same optical design, observatory, and observing 
strategy as BICEP, but will be upgraded to an advanced detector-array 
of superconducting transition-edge sensors (TES). The TES [M] replaces 
the NTD-Ge semiconductor thermistor bonded to each absorbing grid in 
the current PSB design with a superconductor operated near its normal-to- 
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superconducting transition temperature. CMB photons heat up the super- 
conductor, causing its resistance to change enormously, which makes it an 
ideal sensor. 

BICEP requires three electroformed feedhorn antennas per pixel both to 
receive power from the sky and to couple millimeter wave radiation to the 
bolometers, see Figure 11.51 The electroforming process is costly and time- 
consuming, and the feedhorn's dimensions fundamentally limit the pack- 
ing efficiency of the focal plane. The TES methodology has a significant 
logistical advantage: the superconductor and associated components are 
all fabricated photolithographically, resulting in robust, reproducible, mass- 
produced arrays that can subsequently "tile" BICEP-II's focal plane. The 
tiling with planar arrays is extremely efficient — more than five times as 
many TES detectors can be placed in the same focal plane BICEP's 
current NTD-Ge semiconductor array. 

BICEP-II will use TES bolometers with integrated planar antennae |65j 
developed at NASA/JPL and read-out by a time-domain SQUID multi- 
plexer developed at NIST [66]. The "ultimate" large angular scale B-mode 
experiment will be a small array of BICEP-II receivers at the South Pole. 
This array would use current technology to probe the GWB if inflation oc- 
curred at, or slightly below, the GUT-scale, thereby testing all potentially 
observable models of inflation. All of this can be accomplished now at the 
South Pole, for approximately 1% of the cost of a similarly capable satel- 
lite mission, allowing us to obtain perhaps the most enigmatic image ever 
captured: the birth pangs of the Big Bang! 
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